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Abstract. Neutrino emissivity from the electron i/D 
bremsstrahlung in the liquid layers of the neutron star 
crusts is studied. Nuclear composition of matter in neu- 
tron star crusts is considered for various scenarios of neu- 
tron star evolution. For the deep layers of the crust, the 
compositions of cold catalyzed matter, accreted matter 
and hot matter (T > 5 x 10^ K) are shown to be very 
different, and this implies differences in the neutrino emis- 
sivity at given density and temperature. Neutrino-pair 
bremsstrahlung, due to collisions of relativistic degenerate 
electrons in a Coulomb liquid of atomic nuclei, is consid- 
ered. The neutrino energy loss rate is expressed in a simple 
form through a Coulomb logarithm L - a slowly varying 
function of density, temperature, and nuclear composition. 
Non-Born corrections, thermal width of electron Fermi 
sphere, and finite sizes of nuclei are taken into account. 
Implications for cooling of neutron stars are discussed. 



^ ' 1. Introduction 



A neutron star crust extends from the stellar surface to the 
internal core, up to densities of about 10^"* g cm~'^. The 
properties of crustal matter are most important for study- 
ing neutron star physics and observational data (surface 
thermal X-ray emission, X-ray burst phenomena, evolu- 
tion of magnetic fields, etc.). This article has two goals. 
First, we analyze nuclear composition of matter in neu- 
tron star crusts. Second, we reconsider the neutrino-pair 
bremsstrahlung (NPB) in liquid phase of the matter due 
to scattering of electrons on atomic nuclei, 



+ {Z,A)^e + {Z, A) + v + v. 



(1) 



NPB is known to be one of the most powerful mechanisms 
of neutrino emission in matter with density lower than the 
nuclear density, po = 2.7 x lO^** g cm~^. 
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NPB was considered in a number of works (see Itoh et 
al. 1989, and references therein). Let us mention detailed 
studies of Festa & Ruderman (1969), Flowers (1973), Di- 
cus et al. (1976), Soyeur & Brown (1979), Itoh & Kohyama 
(1983), and Itoh et al. (1984a) who analyzed NPB due to 
electron scattering in liquid and crystalline dense mat- 
ter. Pethick & Thorsson (1994) have shown that the band 
structure effects of electrons in solid matter drastically 
suppress the static lattice contribution to NPB. Finally, 
the reconsideration of the neutrino energy loss rate due to 
the electron-phonon scattering in the crystal phase was re- 
cently performed by two of us (see Yakovlev & Kaminker 
1996). 

In the present paper, we will reconsider the problem of 
the calculation of the neutrino losses from NPB in the hot 
neutron star crusts. We will restrict ourselves to the case 
of a liquid phase of the crust (temperature above melting 
temperature at given density). The actual nuclear compo- 
sition of the crust depends on the formation scenario. In 
contrast to previous studies, our calculations will be per- 
formed for realistic models of nuclear composition, which 
are based on specific assumptions about the formation of 
the neutron star envelope. 

The paper is organized as follows. In Sect. 2 we ana- 
lyze physical conditions and nuclear composition of neu- 
tron star crusts. We discuss those differences between var- 
ious neutron star crusts, which are important for the NPB 
emissivities. In Sect. 3 we recalculate the NPB energy 
generation rate for degenerate relativistic electrons and 
Coulomb liquid of atomic nuclei. Contrary to the previ- 
ous works, we include the non-Born corrections and the 
effects associated with thermal width of the electron Fermi 
sphere. We express our results in a simple form - through 
a Coulomb logarithm, which is a slowly varying function of 
density, temperature, and nuclear composition of matter. 
In Sect. 4 we discuss NPB for various models of matter in 
neutron star crusts, and indicate possible implications of 
our results. General formalism of NPB is presented in Ap- 
pendix A, mathematical aspects of the thermal broaden- 
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ing of electron Fermi sphere are considered in Appendix B. 
Dynamic screening of the electron-nucleus interaction is 
outlined in Appendix C. 

2. Composition of neutron star crusts 

2.1. Physical conditions 

Consider dense matter of a neutron star crust in a den- 
sity range 10^ g cm~^ ^ ^ 10^'' g cm~^. If p < pnd — 
(4 — 6) X 10^^ g cm^'', matter is composed of bare atomic 
nuclei (complete pressure ionization) immersed in an al- 
most ideal gas of relativistic degenerate electrons. This 
phase of matter constitutes the outer crust of a neutron 
star. Above the neutron drip density, pnd, some fraction 
of neutrons is not bound in nuclei, and matter is composed 
of nuclei, immersed in the electron and neutron gases. The 
actual value of Pnd is weakly model dependent. It depends 
also on the scenario of formation of the neutron star crust. 
At p > 10^^ g cm~^, the topology of the nucleon distri- 
bution can be very different from the standard one (e.g., 
neutron gas bubbles in nuclear matter with a large neu- 
tron excess, see Lorenz et al. 1993 and references therein). 
Finally, for p > ph — 1-5 x 10^^ g cm~^, nucleons form a 
single, homogeneous phase. Note that the quoted value of 
Ph obtained recently by Lorenz et al. (1993) is significantly 
lower than the standard one, based on earlier calculations 
(Shapiro & Teukolsky 1983). Matter with pnd < P < Ph 
forms an inner crust. 

The state of degenerate electrons in a neutron star 
crust is described by the electron Fermi-momentum pp or 
relativistic parameter x 



Pf = h{3n'^ne)^^^, x = 



^ « 1.009 ' 

mc 



1/3 



(2) 



where jjLf. is the number of baryons per electron, and pg is 
density in units of 10^ g cm~^. The electron degeneracy 
(Fermi) temperature is 



Tf = To (VTT^ - 1), To = — « 5.930 x 10^ K. (3) 

Kb 

We will consider strongly degenerate electrons, T <C Tp. 

The state of the nuclei (ions) is defined by the Coulomb 
coupling parameter 



^5/3 

0.2254 X, 

Ts 



(4) 



where Ze is the nucleus charge, a = [3/(47rnj)]^/^ is the 
ion-sphere radius, and Tg is temperature in 10* K. At 
high temperatures, when F <C 1, the nuclei constitute a 
Boltzmann gas. At lower T (higher F) the gas gradually 
(without any phase transition) transforms into a Coulomb 
liquid. The liquid solidifies (forms a Coulomb crystal) at 
T = Tjn. For a classical one-component plasma of nuclei, 
Tm « 1.32 X 10^^5/3 (p6//xe)^/^ K Corresponds to F = 172 



(Nagara et al. 1987). For high densities and light nuclei, 

the crystallization may be affected by zero-order quan- 
tum vibrations of nuclei (Mochkovitch & Hansen 1979, 
Chabrier 1993). 

2.2. Models of matter in a crust 

If T < 5 X 10^ K ~ To, one can approximate the composi- 
tion of the crust by that calculated at T = 0. Finite tem- 
perature corrections at T ~ 10^ K result in the presence 
of a very small fraction of free neutrons even at p < pnd, 
with a negligible effect on the bulk composition of matter. 

The compositions of the outer and inner crusts depend 
on the scenario of formation of these outer layers. We con- 
sider three models of matter in a neutron star crust. The 
first model is based on the assumption that matter is in 
its ground state (cold catalyzed matter). The second one is 
the model of accreted matter valid when the neutron star 
accreted a sufficient amount of matter during its life. Both 
models are appropriate for T < 5 x 10^ K. In both cases 
the nuclei are assumed to constitute a one-component 
plasma {A, Z), which behaves either as a Coulomb liquid 
(T > Tm) or as a Coulomb body-centered cubic crystal 
(T < Tm). In addition we consider (Sect. 2.4) the model 
of hot matter vaUd for T > 5 x 10^ K. 

Some uncertainties in calculations of nuclear compo- 
sition come from extrapolation of the laboratory nuclear 
physics to very large neutron excesses and huge pressures 
in neutron star crusts. These uncertainties are more im- 
portant for the inner crust, and they increase with increas- 
ing p. 

2.3. Cold catalyzed and accreted matter, T < 5 x 10^ K 

Let us first consider the models of cold catalyzed and ac- 
creted matter in the outer crust (p < pnd)- 

At finite temperature the first model corresponds to 
a complete thermodynamic equilibrium. This assumption 
is certainly valid at T > 10^° K (neutron star birth); its 
validity is extrapolated to the later stages of neutron star 
cooling. This determines the composition at a given pres- 
sure (density). In our analysis of the outer crust, we use 
the model of Haensel & Pichon (1994) based on the most 
recent experimental data on neutron rich nuclei. 

The model of accreted matter represents the case, in 
which matter is not in the ground state. At typical ac- 
cretion rates onto a neutron star in a close binary sys- 
tem, freshly accreted hydrogen-rich matter burns into he- 
lium. The latter, in turn, burns explosively into ^^Ni, 
which transforms eventually into ^^Fe. The products of 
thermonuclear burning are subsequently compressed un- 
der the weight of accreted matter at relatively low temper- 
ature. For instance, at the accretion rate 10^^° Mq/jy 
typical temperature within the crust is a few times of 
10* K. Accordingly the only further nuclear processes in 
the steadily formed 'new outer crust' are electron captures. 
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which eventually lead to neutron drip at pnd — 6 x 10^^ 
g cm""^ (e.g., Haensel & Zdunik 1990). Note, that at the 
accretion rate 10~^° Mq/jy it takes 10^ yr to replace the 
original outer crust, built of catalyzed (ground state) mat- 
ter, by the new one. 

The compositions of the accreted and cold catalyzed 
outer crusts are vastly different. At p ~ 10^^ g cm~^ the 
values of Z and A in the accreted matter arc less than 
half of those in the ground state. For instance, at the neu- 
tron drip we have Z=18, A=56 for the accreted matter 
(Haensel & Zdunik 1990) and Z=35, A=118 for the cold 
catalyzed matter (Haensel & Pichon 1994). The neutron 
drip density, which separates the outer and inner crusts, is 
rather insensitive to the scenario of the outer crust forma- 
tion (cold catalyzed matter or accreted matter) and/or to 
the model of nuclei present in the matter. Irrespectively of 
the model, one gets pnd = (4 — 6) x 10^^ g cm~^ (Shapiro 
& Teukolsky 1983, Haensel et al. 1989, Haensel & Zdunik 
1990). 

Now consider the models of ground state and accreted 
matter in the inner crust (p > Pnd)- In both cases we 
assume that (for p < 10^^ g cm~^) matter is composed 
of spherical, neutron rich nuclei (A, Z) immersed in the 
electron and neutron gases. Because of relatively low tem- 
perature, the only processes involved in the formation of 
the "inner accreted crust" are electron captures, neutron 
emission and absorption, and pycnonuclear fusion. 

For the ground state matter at p > pnd, we use the 
results of Negele & Vautherin (1973). The model of ac- 
creted crust is taken from Haensel & Zdunik (1990). We 
have constructed a specific model of an accreted crust for 
a neutron star which accreted 6 x 10"'' Mq. The bottom 
of the accreted layer is found at p ~ 1.1 x 10^^ g cm""^. 

Both models are developed for a one component 
plasma (single A, Z nuclide at a given pressure). This im- 
plies the onion structure of the crust, with fixed values of 
A, Z within certain pressure intervals. 

Figures 1 and 2 present the relevant parameters of the 
ground state and accreted crusts versus density. The dis- 
continuities in the values of (A, Z) are due to the shell 
and pairing effects in the binding energies of nuclei. No- 
tice a strong dependence oi Z^ /A (which is important for 
NPB, Sects. 3 and 4) on the crust formation scenario for 
10^ g cm-3 < p < 10^3 g cin-3. The values of Z'^/A for 
the accreted crust are about 2 — 3 times lower than for 
the ground state. The differences of the melting tempera- 
ture are even higher (a factor of '--^ 6 at p lO'^^ g cm^'^. 
Fig. 2). The melting temperature of the ground state mat- 
ter for p > 10^^ g cm~^ is > 5 X 10^ K. At such tempera- 
tures the thermal corrections to the composition of matter 
become important. 

Let us notice that in the case of accreted crust at 
p > 10""^ g cm~^, the ion plasma temperature, Tp = 
huip/kB = (?i/fcB)(47rZ^e^ni/mi)^/^ , where rii is the num- 
ber density of nuclei (ions) and mi is their mass, becomes 
larger than the melting temperature of classical plasma. 



This seems to indicate, that corrections to Tm, resulting 

from quantum (zero-point) vibrations of nuclei, might be- 
come significant. However, a quantitative analysis, based 
on the formulae derived by Chabrier (1993) (who corrected 
formulae obtained by Mochkovitch & Hansen (1979) ) 
shows, that even in the most extreme cases, corresponding 
to minima of Tm at p > 10^^ g cm~^ (see Fig. 2), quantum 
corrections to Tm were about a few percent. In the case of 
the ground state of matter we have always Tp < T^, and 
so the quantum corrections to Tm are even smaller. 
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Fig. 1. Nuclear factor Z'^ /A, important for NPB, and the mass 
fraction of free neutrons, Xn (related to mass fraction of atomic 
nuclei, Xa, by Xa = 1 — Xn), versus density p (in g cm~^), 
for two models of neutron star crust. Solid line: ground state 
of matter; dotted line: accreted crust. 



2.4. Hot matter, T>5xlO^K 

For T > 5 X 10^ K, the thermal effects strongly mod- 
ify the nuclear composition. Both shell and pairing ef- 
fects are washed out from the abundances of nuclei. Also, 
reshuffling of nucleons is no longer blocked efficiently by 
Coulomb barriers, and one can assume complete thermal 
equilibrium. We treat hot matter as a mixture of neutrons, 
protons, electrons, positrons and nuclei. In addition, we 
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Fig. 2. Melting temperature (K), versus density (in 
g cm~^), for the same models as in Fig. 1. 



represent light nuclei by the a particles, and heavy nuclei 
by a single species Z) (e.g., Burrows & Lattimer 1984). 

Wc describe the nuclei in hot dense matter, using the 
model of Lattimer & Swesty (1991), with a specific choice 
of the equilibrium incompressibility of nuclear matter, 
it's = 180 MeV. We assume nuclear equilibrium, as well 
as beta equilibrium of matter. The assumption of nuclear 
equilibrium is justified by high temperature. Beta equi- 
librium is adopted for simplicity: a very rapid cooling of 
matter at highest temperatures can produce some devia- 
tions from beta equilibrium. 



In contrast to the nuclear composition, the value of 

Z"^ I A depends on temperature rather weakly (Fig. 4). Its 
density dependence is smooth, because nuclear pairing and 
shell effects are washed out by temperature. 
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Figure 3 shows the composition of matter for T 

5 X 10^ K , 8 X 10^ K, and 1.2 x 10^" K. One 
a temperature dependence for p < 10^'^ g cm""^. 
other hand, for p > 10^^ g cm^'^, the temperature effects 
are negligible. If T ^ 5 x lO*' K, the thermal effects are 
weak and produce a small admixture of free neutrons at 
p < pND- The fraction of free neutrons for p < /9nd rapidly 
increases with growing temperature, and small fractions 
of a-particles and free protons appear at p < 10^" g cm~^ 
(T = 8 X 10^ K). If T = 1.2 X 10^° K, the nuclei evaporate 
completely for p < 10^ g cm^'^, and the mass fraction con- 
tained in free neutrons is significant at any density under 

discussion. For p > 10^^ g cm~^, the temperature effects 3, NPB energy loss rate 
are negligible. 
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Fig. 3. Mass fractions of various particles versus density at 
different temperatures Tg = T/(10^ K) for the model of hot 
matter, described in the text. 
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Fig. 4. The nuclear factor Z'^ /A, important for NPB, versus 
matter density p, for different temperatures of hot matter. 



3.1. General remarks 

Let us consider NPB (Eq. (|l|)) in the case of relativistic 
degenerate electrons in liquid phase of matter of a neutron 
star crust (T > T^). NPB has been studied by many 
authors (Sect. 1). Thus we omit the details and outline 
the derivation of the neutrino energy generation rate in 
Appendix A. 

At the first stage, we adopt the Born approxima- 
tion and neglect energy exchange between electrons and 
atomic nuclei. We will discuss the non-Born corrections 
in Sect. 3.3, and the dynamic effects of nucleus response 
in Appendix C. We will mostly use the units in which 
c = h = = \. 

Let P = (e, p) be the electron 4-momentum in the 
initial state, P' = (e', p') the electron 4-momentum in the 
final state, and Q — (il,q) be the 4-vector of momentum 
transfer to a nucleus. Since the energy transfer to a nucleus 
is neglected we assume $7=0 throughout Sects. 3 and 4 
(although we present general equations in Appendices A 
and C). Let Ki — (a;i,ki) and K2 — (0^2, be the 4- 
momenta of neutrino and antineutrino, respectively, and 
K = K1+K2 = {oj, k) be the 4-momentum of the neutrino 
pair (to ^ uji + LU2 and k = ki + k2). Energy-momentum 
conservation implies 



P = P' + Q + K. 



(5) 



As shown in Appendix A, the neutrino energy produc- 



tion rate Qbic 
gas is 



(ergs cm ^ s in a relativistic electron 



G^pC^ rii 
12(27r)io 



dk ydp Jdp'\U{q)f 



(6) 



where Gf = 1.436 x 10 ''^ ergs cm'^ is the Fermi weak 
couphng constant, Cl ^ + + 2{C'y + C'a) IS ex- 
pressed in terms of vector and axial vector constants (Ap- 
pendix A) and takes into account generation of electron 
neutrinos (Cy and Ca] weak neutral -I- charged currents) 
and also of muonic and tauonic neutrinos {C'v and C a] 
neutral currents), f/(q) is the Fourier transform of the 
electron-nucleus Coulomb potential, J+ is the spin aver- 
aged squared matrix element in the limit of ultrarelativis- 
tic electrons. Furthermore, 



/ = 



1 + exp 



T 



n -1 



(7) 



is the Fermi-Dirac function for the initial electron, /' = 
fie') is the same function for the final electron, and is 
the electron chemical potential. Integration in (6) is to be 
carried out over the domain where > 0. 

It is convenient to express Qsrem in the form (in the 
ordinary physical units) 



Qi 



2-2 

3.229 X 10^^ pi2 Xa -r TIL ergs s'^ cm'^, (8) 



where Xa is the mass fraction contained in nuclei, and 
P12 is density in the units of 10^^ g cm~'^. Numerical ex- 
pression for QBrcm is obtained for the emission of v^, Vp,, 
Vt, C\ w 1.675, with the Weinberg angle sin^ 9iy ~ 0.23 
(see Appendix A). In Eq. (8) we have introduced the di- 
mensionless quantity L. We will see (Sects. 3.3 and 3.4) 
that L ~ 1 is a slowly varying function of density, tem- 
perature and nucleus parameters which has meaning of a 
Coulomb logarithm. Thus the problem reduces to evaluat- 
ing L. Note that if we replace 3.23 x 10^^ L by 2.18 x 10^'^, 
then Eq. (8) transforms into the well known formula of 
Soyeur & Brown (1979). 

3.2. Screened Coulomb potential 

The squared Fourier transform of the Coulomb potential 
screened by the plasma polarization (which enters Eq. (6)) 
can be written as 



\u{^)\ 



2^2 



(9) 



In this case F{q) is a nuclear formfactor which takes into 
account the proton charge distribution in a nucleus, S{q) 
is the static structure factor of ions that describes the ion 
screening due to ion-ion correlations, and e(g) is the static 
longitudinal dielectric function of the electron gas, which 
accounts for the electron screening. 

It is easy to show that the electron screening is always 
static for the conditions of study. The dielectric function of 
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the degenerate electrons was derived by Jancovici (1962). 
In the ultrarelativistic Umit 



Cm 



where 



6y 



hi 



1 



y 



(10) 



3.3. Coulomb logarithm 

Using Eqs. (6, 8) and we come to the general expression 
for the Coulomb logarithm (here and below again h = c = 
fee = 1): 



L 



189 



2ii7r9T6 
X fil-f')LuJ+ 



dk dp dp' 



1 Siq)\F{q)\'^ 



q^ee' \e{q)\ 



(14) 



y^T, > fcTF=PF 



1/2 



fcxF being the electron static screening momentum, and v-p 
the electron Fermi velocity. In the long-wavelength limit 
{q <C pf), one has 



g2 



which corresponds to the Debye screening. Actually this 
approximation appears to be quite satisfactory although 
we have used exact dielectric function (10) in numerical 
calculations. 

For a spherical nucleus with a uniform proton core, 



Fiq) 



[sin(gro) - qro cos(gro)] . 



Let us simplify the integration. Since the electrons are 

(11) strongly degenerate, the main contribution to L comes 
from those electron transitions in which the electron mo- 
menta p and p' lie in the narrow thermal shell around 
the Fermi surface, |e — /x| < T and \e' — fJ.\ < T. Let 
qs ^ 1/a be the typical screening momentum (Sect. 3.2). 
We assume that qs pf to simplify our analysis. 

Let q = Qt + , where corresponds to purely elastic 

(12) Coulomb scattering while q^. takes into account inelastic- 
ity. Here the inelasticity means that the length of electron 
momentum changes slightly due to the Coulomb interac- 
tion although the energy transfer to the nucleus is absent. 
The process is kinematically allowed since two neutrinos 
are also involved in energy-momentum conservation. Let 
us also introduce the vector p" = p — q^ = p' + q,, + k 
which is directed along q^ but has the same length as p. 

Note that qt = 2psm{6/2) « 2pF sin(6'/2), where 9 is 
^2^3-) an angle between p and p". From geometrical considera- 
tion we obtain 



Here tq is the core radius which can be smaller than the 
nucleus radius, for neutron- rich nuclei. One has F(q) ~ 1, 
for point-like nuclei. 

The structure factor S{q) in (||) has been calculated by 
many authors (e.g., Itoh et al. 1983) for a one-component 
classical plasma of ions with the uniform electron back- 
ground. In a strongly coupled plasma (F > 1), S{q) can- 
not be evaluated analytically. Its general feature is that it 
is greatly suppressed in the long- wavelength limit, qa < 1 
(a is defined in Eq. (|^)). Thus S{q) produces the screening 
of the Coulomb interaction with the screening momentum 
~ 1/a. 

Note that actually the structure factor can be influ- 
enced by response of neutron gas, which surrounds atomic 
nuclei in dense matter, to the motion of nuclei. This effect 
has not been considered so far in the literature. We will 
use the conventional structure factors (Itoh et al. 1983), 
assuming that the inclusion of this effect can be reduced 
to proper choice of the fraction of free neutrons and the 
number of neutrons bound in nuclei (i.e., proper determi- 
nation of the nucleus mass number). 

The electron and ion screenings suppress the Coulomb 
interaction for small q while the nuclear formfactor reduces 
the interaction for large q >, I/'Tq- For typical parameters 
of dense stellar matter, the ion screening is more efficient 
than the electron one (^tf < a~^). 



q'-i<lt+<lr)'-qf+<l'r~1^[^)- (15) 

The neutrino-pair momentum k can be presented as 
k — kt -|- kr, where k^ and kf are the orthogonal vector 
components parallel and perpendicular to p", respectively. 
Strong electron degeneracy implies that q^ and k are much 
smaller than p^ although qt can be comparable with pp for 
large-angle electron scattering events. Therefore, NPB is 
accompanied by nearly elastic Coulomb scattering. Then 
the neutrino-pair energy is 

2 /2 

£ — £ 

ijJ = — « kr + qr- (16) 

£ + £ 

The condition K'^ = cj^ — k^ « fcg — fc| > requires 
ko > kt, where fcg = qr{2ijj — q^), or > 0, w > qr/2. 
Furthermore, we set dp = dfie^de, where dfi is solid angle 
element in the direction of vector p. The integration over 
£ under strong electron degeneracy is standard: 

The integration over p' can be replaced by the integration 
over q, with dq = 2'nqt Aqt dq^- Note also, that dfc^ = duj. 
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This yields 
189 



2ii7r9T6 
dui 



8n' 



dqr / dqt qt 



q,./2 



o/T _ I 



dkt kt 



27r 



where if is an azimuthal angle of k with the polar axis 
along p". The expression for J_|_ is given by Eq. (A14) in 
Appendix A. Using this expression and the relation q^k = 
—2ppkrSm^{6/2) + ppkt sin 9 cos (p, we can integrate over 
ip and kt and present L in the form 



L = 



1 

f 



2pF 



dqt qt 



dqr 



S{q)\F{qW 



use the well known method of soft photons (Berestetskii 
et al. 1982) and claim that the NPB rate is proportional 

S{q)\F{q)\' 

to the product of the cross section of the elastic electron 
(7*|e((7)p - nucleus scattering and the probability of the neutrino- 

pair emission. The pair emission involves weak interaction 

df J+, (18) 

and cannot be affected strongly by the Born approxima- 
tion. Hence the main impact of the Born approximation is 
on the elastic scattering cross section. We expect that the 
non-Born correction factor is i?NB(9t) = cr(gt)/(7Born(9t), 
where a{qt) is the exact elastic cross section with the mo- 
mentum transfer qt, and crBorn('?t) is the Born cross sec- 
tion. The factor i?NB('?t) has been calculated by many 
authors (e.g., Doggett & Spencer, 1956), and its inclusion 
into the Coulomb logarithm is straightforward. 



X RT{y,u) RisiBiqt), 



(19) 3.4- Low- and moderate-temperature cases 



where y ~ qt/{2pp) — sin(6'/2), u = qr/T, and q is given 
by (^5|). Here we introduce the factor i?NB to take into 
account deviations from the Born approximation (see be- 
low) . The integration over describes the thermal effects, 
i.e., a weak inelasticity in Coulomb scattering due to a fi- 
nite width of the thermal shell around the electron Fermi 
surface. In Eq. (19) these effects are controlled by the func- 
tion 



Rt{v,u) 



,ui{ir) „2 ^ 

dv \v — — 

e^ - 1 V 2 




(20) 



where v — lo/T. The thermal effects introduce an ad- 
ditional screening of the Coulomb interaction with the 
screening momentum T. The properties of Rriu, u) are 
analyzed in Appendix B, where we present also a conve- 
nient fitting expression for this function. 

Finally, let us discuss the non-Born correction i?NB 
in (19). Since dense stellar matter contains heavy nuclei 
(Sect. 2), the Born approximation used in all previous 
works is not very accurate. Exact calculation of the NPB 
rate beyond the Born approximation is difficult but we 
propose an approximate treatment of the non-Born terms, 
which is explained below. 

A process of NPB consists of two stages: electron scat- 
tering on a nucleus, and neutrino-pair emission. As seen 
from the above results the electron scattering is much 
stronger (i.e., it is accompanied by larger momentum 
transfers) than the neutrino-pair emission except possibly 
for the small-angle scattering. However one cannot expect 
large deviations from the Born approximation for small- 
angle scattering (Berestetskii et al. 1982). Thus we can 



We have reduced the problem of calculating the Coulomb 
logarithm to a two dimensional integration in (19). One 
can distinguish two cases: the low-temperature case when 
temperature is lower than the Coulomb screening momen- 
tum, T <C and the moderate-temperature case when 

All previous studies of NPB from the relativistic degen- 
erate electrons have used the approximations appropriate 
to the low-temperature case. In this case the main con- 
tribution into L comes from the values of qt ^ T. Then 
Rt is a sharp function of which decreases rapidly with 
increasing qr ^ T (Appendix B). Therefore, we can set 
qr = in the remaining functions under the integral (19), 
and the integration over qr is performed with the aid of 
Eq. (Bl): 



L 







qMqW 



where q — qt, y ^ q/ (2pf) and 
2y' 



Rciv) - 1 



1 



y 



■ln{y). 



(21) 



(22) 



Equations (8) and (^l|) with i?NB — 1 reproduce the famil- 
iar NPB rate for degenerate and relativistic electrons. For 
example, these equations can be obtained from the results 
of Itoh & Kohyama (1983), taking into account that 



yt. 



yx + ^ 



yx-S. 



y 



(23) 



where ^ = I — x'^ {I — y^) . The quantity -Fiiquid intro- 
duced by Itoh & Kohyama (1983) is equal to (2/3)i, and 
their function hiq) = (8/3) R^{y). 

We see that L has, indeed, the meaning of the Coulomb 
logarithm. If the screening of the Coulomb potential were 
weak {S{q) — e{q) = 1) and i?c ~ 1 (which is true 
for q <C pf), then L would acquire a familiar logarith- 
mic divergency at small q. The divergency is eliminated 



8 



P. Haensel et al: Neutrino pair bremsstrahlung 



due to the screening functions (Sect. 3.2). The function 
Rc, Eq. (p^, introduces additional screening which comes 
from the squared matrix element of the NPB reaction. 
We have i?c — ^ 1 for small-angle scattering (y ^ 0) , while 

w {l — y^)/2 for backscattering {y — > 1). Therefore 
Rc suppresses the backscattering (just as F{q)) which is 
natural for relativistic electrons (Berestetskii et al. 1982). 

The low-temperature Coulomb logarithm L is ex- 
pressed as a simple one dimensional integral (pl|), which 
is easily computed (Sect. 3.5) for any parameters of dense 
stellar matter. 

In the moderate-temperature case {Qs 1^ T ^ Tp), one 
should deal with two dimensional integration (19) which 
is also easy once the thermal function i?T(y,u) is known 
(Appendix B). 

3.5. Numerical results 

In the degenerate electron gas and Coulomb liquid of 
atomic nuclei, the Coulomb logarithm depends actually on 
four parameters: on nuclear charge number Z that defines 
the ion screening (Sect. 3.2), on the proton core nuclear 
radius ro (see ([T^)), on density p and T. The ion coupling 
parameter F, Eq. (^, is expressed through p and T. 
In the ultra-relativistic limit (x 3> 1) the four parameters 
can be replaced by three dimensionless parameters 



L{Z,f],t), 7] 



a 



t 



kBT 
2pFC 



T 

27>' 



(24) 



where a is the ion-sphere radius (see Eq. (Q)), and Tp is 
given by Eq. (||). 

We have calculated L from Eq. (19) for Z = 10, 20, 
30, 40, 50, and for wide ranges of ry < 0.2 and t < 0.1 
in the Coulomb liquid (T > T„„ i.e., i > 1.1 x IQ-'^Z^^^, 
Sect. 2.1). The structure factor S{q) has been taken from 
Itoh et al. (1983). Typical results are presented in Figs. 5 
and 6. 

At small t, our values of L tend to those obtained in 
the low-temperature approximation (cf., curves 3 and 4 
in Figs. 5 and 6). At higher t, the deviations from the 
latter approximation are quite pronounced. In the low- 
temperature approximation, L depends on T only through 
the structure factor S{q) (which slightly varies with F), 
and L grows slowly with t due to weakening of the ion 
screening. An inclusion of the thermal effects (Sects. 3.3 
and 3.4) leads to a stronger and non-monotonic temper- 
ature dependence of L which can be explained as fol- 
lows. The thermal effects are described by the function Rt 
which, generally, introduces an additional 'thermal screen- 
ing' of the Coulomb interaction and tends to suppress L 
and the NPB rate. However the thermal effects act within 
the thermal shell near the electron Fermi surface. Devi- 
ations of the electron momenta from the Fermi sphere 
can decrease the momentum transfers q, Eq. ([l5|), for a 
given qt. Since the NPB rate, Eqs. (8), (19), involves q^"^ 
(squared Coulomb potential), the decrease of q causes a 



noticeable growth of L with t at low t (Figs. 5 and 6). 
With increasing t, the 'thermal screening' itself becomes 
more important and suppresses L at higher t. 

Figures 5 and 6 also show variation of the Coulomb log- 
arithm with Z . The Coulomb logarithm increases with Z 
which is evidently the effect of the ion screening (qs/pF ~ 
Z^^^^): the smaller the ratio qs/pF, the larger the reaction 
rate. Similar effects in the low-temperature case are well 
known (e.g., Itoh & Kohyama 1983). 



1.4 
1.3 
1.2 
1.1 
1.0 
0.9 
0.8 
0.7 



[ Z=30 ^ 


" ^ 1 ' 




" 2 i 

3 ; 




^ 4 : 







0.02 0.04 0.06 0.08 
t 



0.1 



Fig. 5. The NPB Coulomb logarithm L vs. dimensionless tem- 
perature t (^) for Z = 30 nuclei. Curve 1 - point-like nuclei 
(77 = Vc/a = 0); 2 - finite-size nuclei with rj — 0.2; 3 - same as 1 
but in the Born approximation; 4 - point-like nuclei, non-Born 
corrections and thermal eff'ects are neglected ('standard ap- 
proach' of previous works). 



The calculations show that finite size of nuclei (nu- 
clear formfactor) becomes significant for 77 > (0.1 — 0.2). 
With increasing 77, the Coulomb logarithm noticeably de- 
creases since the nuclear formfactor introduces an effective 
screening of the Coulomb interaction (Sect. 3.2). Finite- 
size effects are negligible in the outer crust of a neutron 
star, where 77 ^ 1, but they are strong in the inner crust, 
where the atomic nuclei occupy a substantial fraction of 
volume. 

Finally, let us emphasize the importance of the non- 
Born corrections. According to Doggett & Spencer (1956), 
the non-Born correction factor i?NB increases the electron 
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t 



Fig. 6. Same as in Fig. 5, but for Z=50. 

- nucleus scattering cross section, i?NB > 1: and the in- 
crease is larger for higher Z. Accordingly, an inclusion of 
this factor amplifies quite noticeably the Coulomb loga- 
rithm as seen from Figs. 5 and 6. For Z > 30, this effect 
is even more important than the thermal effect described 
above. 

Comparing our improved values of L which include the 
thermal and non-Born corrections (curves 1 in Figs. 5 and 
6) with the non-corrected values (curves 4), relevant to 
the "traditional" results, we conclude that our corrections 
amplify L and the NPB rate, typically, by a factor of 1.5 

- 2.5. 

For a practical use, we can propose an analytic fit: 

L ^ A/B^/\ (25) 
A = 0.269 20t + 0.0168Z 0.00121r? - 0.0356Zr? 

+ O.OmzH + l.54Ztr], 
B = 1 + 180t^ + 0A83tZ + 2QtZif + 4.31 x IQ-^Z^. 

This formula reproduces all calculated values oi L {Z < 
50, t < 0.1, T] < 0.2) with the mean error of about 1%, 
and with the maximum error of 2.3% at Z = 10, t = 0.001 
and rj = 0.2. 

Thus the NPB energy loss rate QBrem can be easily 
calculated from Eqs. (8) and (25) for any model of dense 
matter (provided the values of p, T, rie, Z, A, Xa, and 
To are specified). If matter consists of nuclei of various 
species {A, Z), one should sum over the species in Eq. (8). 



Actually, NPB from a multi-component mixture of nuclei 

deserves a separate study. We expect, however, that our 
results (obtained for one component plasma of nuclei) can 
be used, at least semi-quantitatively, for the mixtures as 
well. 

Note that Cazzola et al. (1971) and Munakata et al. 
(1987) considered the NPB of non-degenerate and weakly 
degenerate electrons. In this case the thermal shell washes 
out the Fermi surface, and the thermal effects are natu- 
rally implanted in the equations. However the above au- 
thors analyzed ideal plasma of ions F ^ 1 (neglecting 
the ion screening) and weak degeneracy (T/Tp > 0.3) - 
the conditions which are not relevant for applications and 
opposite to those studied in our work. 

Figure 7 shows the density dependence of the NPB 
energy loss rate .BBrem = QBrem/p (logarithmic scale) for 
three models of dense matter (Sect. 2) at several T. For 
the ground state and accreted matter, the NPB is pro- 
duced by electron scattering on nuclei {A, Z) of one species 
(Sect. 2). In the case of the hot matter, we take into ac- 
count the contribution of nuclei, protons and a particles. 
Note that self-consistent models of accreted matter corre- 
spond to T > 10^ K (e.g., Miralda-Escude et al. 1990). We 
use such a model at higher T for illustrative purpose, to 
show how possible scatter in nuclear composition affects 
the NPB rate. The rates for accreted and groimd state 
matters are shown at T < 5 x 10^ K, while the rate in 
the hot matter is displayed at T > 5 x 10^ K (Sect. 2). 
If T ^ 5 X 10® K, all three models arc possibly not very 
accurate. We present three curves at T = 5 x 10® K to vi- 
sualize the effects of nuclear composition. The NPB rates 
at this temperature are shown separately in Fig. 8 in nat- 
ural scale. It is seen that they differ by a factor of 2 - 
4. 

At T = 10® K and 2 x 10® K, in the considered density 
region 10® — 10^^ g cm ~^ the matter may be in liquid and 
crystal phases. The quantity .BBrem in the crystal phase 
was calculated using simple approximations obtained re- 
cently by two of the authors (Yakovlev & Kaminker 1996) 
for neutrino pair emission due to the electron-phonon 
scattering. The ground state matter is crystallized for 
p > 2.35 X 10^° g cm -3 at T = 10® K (not for all 
densities, see below) and for p > 1.3 x 10^^ g cm ~^ at 
T = 2 X 10® K. The accreted matter is crystallized for 
p > 2.14 X 10^^ g cm (not for all densities) only at 
T = 10® K. The melting temperature for the accreted 
matter is smaller (Fig. 2) due to lower values of Z, and the 
crystallization occurs at higher p (if T is fixed). It should 
be noticed that there are two types of jumps of the NPB 
rate for the ground state and accreted models. The jumps 
of the first type result from the jumps of Z and A (Sect. 2, 
Figs. 1, 2). Such jumps take place at T = 10® K, 2 x 10® K 
and 5 X 10® K (Figs. 7, 8). The other more pronounced 
jumps come from transitions from liquid phase to crys- 
tal and back at some densities. One may see single 'back' 
transition from crystal to liquid for the ground state mat- 



10 



p. Haensel et al: Neutrino pair bremsstrahlung 



ter and a few ones for the accreted matter at T = 10^ K 
(Fig. 7). The 'back' transitions come from nonmonotonic 
character of melting temperature due to the jumps of Z 
(see Fig. 2). When T is fixed and p > 3 x 10^ g cm~^, the 
NPB rate for the ground state matter in the hquid phase 
is, generally, several times higher than for the accreted 
matter due to larger values of Z (Sect. 2). 

Note, however, the striking general similarity of the 
NPB rates in the liquid and solid phases (Fig. 7). It re- 
sults from the similarity of the expressions for the NPB 
rates (cf. Eq. (21) with Eqs. (13) and (14) of Yakovlev 
& Kaminker 1996) and reflects common properties of 
strongly coupled Coulomb liquid and high-temperature 
Coulomb crystal. 

The NPB rate in the hot matter for T = 1.2 x 10^° K 
(Fig. 7) is broken at densities p < 2 x 10* g cm^'^ at which 
the electron degeneracy becomes low and our results are 
invalid. Note a sharp drop of the NPB rate at p < 10^ 
g cm^'^ for this T. It occurs due to dissociation of nu- 
clei (Fig. 3): the contribution of nuclei into the NPB be- 
comes negligible but a smaller contribution of protons and 
a-particles is available. The drop clearly indicates the im- 
portance of coherence effect (QBrcm oc Z"^) which amplifies 
significantly the NBP rate in the presence of high-Z nu- 
clei. For other T and p, the contribution of protons and 
a-particles is insignificant. 

If Z, A, Xa, and L were independent of density, the 
NPB rate i?Brcm = QBicm/p would also be density inde- 
pendent (at fixed T). However, as seen from Figs. 7 and 
8, the NPB energy loss rate i?Brcm decreases mainly with 
p at p > 10^" g cm"-^. This decrease is explained mostly 
by the lowering of Z'^/A, (Figs. 1 and 4), fraction of nuclei 
Xa, and of the Coulomb logarithm L. 

4. Conclusions 

We have analyzed nuclear composition of neutron star 
crusts (Sect. 2) for three models of dense matter: cold 
catalyzed matter, accreted matter, and hot matter in ther- 
modynamic equilibrium. The first and second models are 
valid for T < 5 x 10^ K, while the third one is appropri- 
ate for higher temperatures. Nuclear composition {A,Z), 
melting temperature and other properties of dense matter 
are different for these models. 

In Sect. 3 we have reconsidered the neutrino-pair 
bremsstrahlung (^ of relativistic degenerate electrons due 
to scattering on atomic nuclei at T > Tm. We have shown 
that the neutrino energy generation rate can be expressed 
(Eq. (8)) through a Coulomb logarithm L, which varies 
slowly in stellar matter. We have obtained a general ex- 
pression (19) for L, which takes into account two effects 
neglected in previous works. First, it includes the non- 
Born corrections, and, second, it describes the thermal 
effects associated with the finite width of the shell around 
the electron Fermi surface. The thermal effects are shown 
to be quite important at moderate electron degeneracy. 




9 10 11 12 13 

logio/? (g cm ) 

Fig. 7. The NPB energy loss rate iJercm = Qsrcm/p (in 
erg s~^, logarithmic scale) vs. density (in g cm~^) for 
Tg=l, 2, 5, 8, 12 (figures near curves). Solid line - hot mat- 
ter; long dashes - ground state matter, short dashes - accreted 
matter. 

and they lead to the appearance of the new 'moderate tem- 
perature regime' (Sect. 3.4) overlooked in previous stud- 
ies of NPB from relativistic degenerate electrons. We have 
calculated the Coulomb logarithm numerically for possi- 
ble parameters of dense matter in neutron star crusts, and 
found a simple analytic fit (25). 

Note that NPB of electrons in a Coulomb liquid of 
atomic nuclei is similar to NPB at T < produced due 
to Coulomb scattering of electrons on 'charged impurities' 
- nuclei {Ai^p, Zi^p) immersed accidentally in a lattice of 
bulk nuclei {A,Z). When the nuclei crystalhze (T < Tm) 
but temperature does not drop much below Tm, the most 
important is the neutrino pair generation due to electron 
phonon scattering in solid matter (this follows from the 
recent results of Pethick & Thorsson 1994). The neutrino 
energy generation rate drops rapidly with decreasing T at 
T < T,„ (e.g., Itoh et al. 1984b, Yakovlev & Kaminker 
1996), and the Coulomb impurity scattering can be dom- 
inant even for not very low T. The energy generation rate 
Qimp is obviously given by the expression similar to (8) but 
with Hi replaced by the impurity number density riimp, and 
with Z^ replaced by {Zi„^p — Z)^. The expression for Qimp 
contains a Coulomb logarithm iimp which should be quite 
similar to the Coulomb logarithm (^ij) in liquid matter 
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Fig. 8. Same as in Fig. 7, but for T = 5 x 10^ K (linear scale). 
The plotted quantity is ^Brem.io (the loss rate ^Brem in units 
of 10^° ergs s"^ g"^). 



(although the thermal effects become unimportant). The 

major difference is in the Coulomb screening. In solid mat- 
ter, the screening is described by the impurity structure 
factor Simpiq) which takes into account impurity distri- 
bution over lattice sites. If the impurity correlation length 
is large, the main contribution into the screening comes 
from the electrons. 

The results of the present article can be useful for nu- 
merical modeling various phenomena related to the ther- 
mal evolution of neutron stars. First of all, we should 
mention cooling of young neutron stars (of age t < (1 
10^) yrs) where internal thermal relaxation is not achieved 
(Lattimer et al. 1994). The thermal relaxation is accom- 
panied by the propagation of the cooling wave from the 
interior to the surface. The associated variations of surface 
temperature are, in principle, observable. The dynamics of 
thermal relaxation is sensitive to the properties of matter 
in the neutron star crust, particularly to nuclear composi- 
tion and neutrino generation mechanisms. In addition, the 
above results arc useful for studying thermal evolution of 
accreting neutron stars (see, e.g., Miralda-Escude et al. 
1990, and references therein). 
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Appendix A: NPB matrix elements and energy 

loss rate 

Using the notations introduced in Sect. 3.1, the NPB energy 
loss rate in the Born approximation can be generally written 

as (/i = /CB = c = 1): 



where 



dp / dp' / dk,y / dkj, 



(27r)" 

X 5is-e' -uj)cjfil-f')W, 



2 {2oj^){2ujl){2e){2s') 



E 



(Al) 



(A2) 



Afp is the squared matrix element. Summation is over the 
electron spin states a before and after scattering and over neu- 
trino flavors (fe, v^, Vr). When the neutrino energies are much 
lower than the intermediate boson mass (~80 GeV) the stan- 
dard approach yields 



Ei^i' = / 



dq|f/(q)|^5(p-q- 



k) 



X Tv{KiO'^ K20'^) 

X Tr [(P' + m)Lc{P + m)L^] , 

0« = 7"(l + 7'), 

= r,G(P - Q)7° + 7°G(P' + Q)ro 
P + m 



G(P) = 



(A3) 

(A4) 
(A5) 



Here f/(q) is the Fourier transform of the Coulomb potential 
(Sect. 3.2), m is the electron mass, G(P) is the free-electron 
propagator, 7" is a Dirac matrix, upper bar denotes Dirac con- 
jugate, and P = PQi7" (Berestetskii et al. 1982). Furthermore, 
Cv and Ca are, respectively, the vector and axial vector weak 
interaction constants. For the emission of electron neutrino 
(charged -|- neutral currents), one has Cv = 2sin^0w 4-0.5 
and Ca = 0.5, while for the emission of muonic or tauonic 
neutrinos (neutral currents only), C'y — 2 3111^ 6w — 0.5 and 
C'a = —0.5. Here Ow is the Weinberg angle, sin^ Ow — 0.23. 
Using the well known identity 



dk. JdKs^^Hi 



we obtain 

^Brem — 



(K -Ki 
+ 2K'^K^ 



K2) 



(A6) 



12(27r) 



^Idkldp Jdp'McOf 



J = Y^{K°'K'^ - K^g°''') 

X Ty[{P' + m)L^{P + m)Lp\ 



(A7) 



(A8) 
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The integration in (A7) is to be carried out over the domain 
where > 0; g"'^ is the metric tensor. 
Let us introduce the notations 



u = [P' + K f ~rn ^ 2P'K + K'^ , 
V = -{P - Kf ^2PK - K'^. 

Tedious but straightforward calculations yield: 

J = cIj+ + c'^_j-, 



(A9) 



i), 



cl = J2^cl - cl), 



where 

\v u) 

+ fi_l) \k'UpP' -8ee' + . 

+ K' - 4Lj(e - e')) - "^m'to'] 

+ ^K'\7n'' - K'')Uee' -2PP' 

+ 2m^ +4e(x) - K') 

+ ^K'\77i' - K')Uee' -2PP' 

+ 2rn -Ae'uj-K') 



+ A:K (ra - _K ) — ■ 



(AlO) 



+ A r4A''^(pp' - m^)(2ee' - PP' - ni') 
uv 

- 277i''K'{2e£' - PP' + rn') 

+ 4uj'{77i' + K')PP' - 2m^a;^(2m^ + K') 

+ {3e' +3e'^ +2ee' -2PP' -77i' -Lu')] (All) 



and 



Vi; 11/ 

\v uJ 



_ ^7n'K'{4££' -2PP' 

+ 2rn + 4etj - 7^^) 

_ i|m27f2(4ee'-2PP' 

+ 2m^ -4e'aj- A"^) 



12m^74:^ 



8m' 



+ ri^2(^2ee' + 2PP' - 27n' + 7^^) 

uv ^ 

- 2lj^(2PP' - 2m^ + . 



(A12) 



Equations (A7), (AS) (All) and (A12) determine the NPB 
rate for any degree of electron degeneracy and relativism. The 
equations yield the solution of the problem in a compact and 
convenient form which allows one to consider easily different 
limiting cases. In a hot, non-degenerate plasma, where a sig- 
nificant number of positrons can be present in addition to the 



electrons, the contribution from positrons should be added to 
Qfirem iu the Straightforward manner. 

Analyzing Eqs. (All) and (A12) one can show that the 
main contribution into NPB of relativistic and degenerate elec- 
trons comes from J+. Accordingly, we keep only this term in 
Eq. (8). 

In our case, the neutrino-pair momentum is thermal, 
k < T <^ Pf but the elastic electron - nucleus momentum 
transfer can be much larger, qt ^ pf- First consider the case 
when g ^ T or g 3> fc. Then we can put q ~ p — p', and the 
Coulomb scattering is, to a very good approximation, elastic. 
A careful analysis of various terms in (All) shows that the 
main contribution comes from the term 



J+ = 



8K'{PP')[2ee' - jPP')] 
{P'K)(PK) 

{e'uj — p'k)(ecij — pk) 



(A13) 



Here we have adopted the natural approximation: u — 2P'K, 
V = 2PK. In addition, we can set e = e' = /i, p = p' = 
Pf and (qp)/p = (f /2p (see Sect. 3.3), in smooth functions 
of electron energy and momentum. One can easily show that 
(££')^-(PP')^ [l-aV(V)], u « 2pqr and t; « 2p(g, - 

(qtk)/p), where we use the notations of Sect. 3.3. Finally, for 
the case of T <C q we obtain 



J+ — 8 — (lu 

qr 



k') 



4p2 



kqi 

P 



(All) 



Now consider J+ in the approximation of small-angle scat- 
tering [q <C Pf). Then k can be comparable to q, and more 
terms should be kept in J+. An analysis of Eq. (All) shows 
that we must retain the terms 



-[(ee ) 



tit; 



(PP 



4K' 



- + - 

\U V J 



+ 167^2 fi-i) (ee' + pp') 

~ 8K\ee' +ii>Y>'){\ + \ 
16K^ 



+ 



-(£ + £ 



(A15) 



Using energy-momentum conservation and the inequalities 
k <C p and q <^ p and introducing the notations of Sect. 3.3, 
we obtain 



LO q,- ~\- K, 
M ~ u ~ 2pqr 



K' 

J+ 



ki 



LU — 

2 

2 i2\<lt 



kf) 



k' 



(A16) 



Let us compare Eqs. (All) and (A16) for J+ obtained, re- 
spectively, in the domains of large and small q. Since both 
domains overlap, we ca n des cribe J+ accurately for all q. How- 
ever, we see that Eq. (A14) coincides formally with (A16) in 
the domain of small q, where the terms q' / (4p^) and kqt/p are 
negligible second-order corrections. Thus we use Eq. (A14) for 
all q in Sect. 3.3. 
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Appendix B. Thermal function RT{y,u) 

According to Sect. 3.3, the effects of inelasticity in the electron- 
nucleus scattering are described by the thermal function 
RT{y,u) given by Eq. (20). 

Integrating over it we obtain 



Appendix C: Dynamic screening 

In principle, the ion screening of the Coulomb interaction can 
be dynamic. For a strongly coupled plasma of ions (F > 1), 
the dynamic screening can be introduced through the dynamic 
structure factor S{q, fi) which satisfies the relationship 



du Rriy, It) = 1 



■\ny 



l-y2 
1 y < 1, 



In the limits of small and large u we have 



RT{y,u) 



63 1-2/2 



for y ^ u, 



S{q) = 



f.(.n). 



(CI) 



Here, hO, is the energy transferred from electrons to nuclei. 

(Bl) 

Typical frequencies f2 of the dynamic response are Q, ^ ujp, 
where ujp = ^iiZ'^e^ni/mi is the ion plasma frequency, and 
rrii is the ion mass. One can see that TiUp/kB < Tm <C If in 
dense stellar matter under consideration. In order to include 
the dynamic screening, we introduce Q into the energy conserv- 
(B2) ing delta function in Eg. (Al), and rewrite the static structure 
factor S(q) in the form (CI). Then the Coulomb logarithm for 
the relativistic degenerate electrons reads (cf. Eq. (14)) 



RT(y,u) 



189 



327r6 

If t/ < 1, Eq. (20) yields 
189 



u^e for w (1 — j/2) S> 2 j/'^ 



RT{y,u) ^ Rt(u) 



,. , d« 



(B3) L = 



189 



21l7j-9r6 
X /(1-/V^4 



dt7 

2-K 



dk dp dp 



, g(g,n)|F(g)|= 
ee'g4|e(g)|2 



(C2) 



This equation can be simplified as in Sect. 3.3. Let us introduce 
the same quantities q = qt + q, , k = kt + k,. and p" = 
p— qt = p'+k+qr , where |p"| = |p|. Using energy-momentum 
conservation and the inequalities k <^ p, qr <^ p we obtain 



The asymptotic forms are: 



= ^r(5)C(5) l»u»y^ 



Rt{u) 



189 
327r6 



2 —u/2 ^ , 



(B5) 



where (^{x) is the zeta function. Finally, for (1 — y) <C 1 we 
have 



0, + UJ ~ gr + kr 

u = 2P'K + ^ 2p{lj - kr) = 2p(g^ - fi), 

V = 2PK ~ ^ 2p{qr -Q- qtk/p), 

r^2 2 72 i2_72 j2 

K — U) — fc^ — fej ~ Kg — fct . 



(C3) 



Here ko — {qr — Q,){2ijj — qr + fl), and the condition > 
yields 2uj > qr — Q > 0. li q ~ qt and q is much larger than qr, 
Q, and uj, we have (cf. Eq. (A14)) 



RT{y,u) 



189 
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{l-y)u 



/2 



av (v I. 

e" - 1 ^ 2' 



J4 



8(fcg - fc?)g? 



(g^ - r2)(gr - n - qtk/p) 



4p2 



(C4) 



(B6) On the other hand, using the approximation of small-angle 
scattering (g <^ p) and (A16), we obtain (cf. Eq. (A15)) 



We have calculated RT{y,u) numerically for wide ranges 
of y and u (10 points of y from y=0 to 0.9; 10 points of it from 
it=0.1 to 30). We have found the following fitting expression 

RTiy,n) = u{l~y') g, 

F = 0.3058 + 4.3311/2 -hFlue''''^ 

G = 1 + 58.051/2 +Giiie", 

Fi = 2.949 - 2.9631/2 + uF2, 

F2 = 0.7184 - 0.85651/2 ^ o.06019ii 

+ 0.07671^2^ _ 0.00077711*2 1/, 

Gi = 9.797 - 6.502y2. 



(B7) 



The mean error of the fit is 3.6%, and the maximum error of 
8.1% takes place at y=0.5 and it=4. This expression has been 
used in calculations presented in Sect. 3.5. 



J4 



8(fco 



(gr - 0)2 



[gi2 + 2^l{qr - n)] 



(C5) 



Combining Eqs. (C4) and (C5), we can propose the following 
interpolation 



J+ 



8(fcg - fc2) 



(g^ - Q){qr ~ Q ~ qtk/p) 
X [g2 + 2n{qr - ft)] . 



4p2 



(C6) 



Calculations similar to those in Sect. 3.3 allow us to express L 
in the form analogous to (19). Using (C6), we come to a simple 
equation for the Coulomb logarithm: 



L = 



27 



2pp 



dqt qt 



^ [g? + 2Q{qr - m 



X ^M^i?.(y,-,-)i^-B(gO. 



(C7) 
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Here q = \J ^ gP — q^iqr/p), and we introduce a new func- 
tion 



RT(y,u, w) 



189 
167r« 

1 



'1 

du 

■"0 

V — vo 



v{v + w){v — Uo)^ 



3i;i(l-y2) 



dv 



v(v + w) 



v-vo- ^vi{l~y^) 



(C8) 



with y = qt/{2pF), u = qr/T, w = Q/T, v = lo/T, vq = 
{u — w)/2 and Vi = vo/y^- This function describes the thermal 
effects (Sect. 3.3) in the presence of the dynamic screening. 
Eq. (C8) yields 



du RT{y, w + u ,w) 



My) 



63 

87r6 



dv 



«^(u + w) 



(C9) 



where Rc(y) is defined by (^2|). Consider the low-temperature 
case {T <^ qs) and take into account that RrijJ, u, w) decreases 
exponentially with increasing (g^ — ^) at (g^ — Q,) > T (see 
Sect. 3.4). Then Eq. (C7) reduces to 



63 
87re 



2PF 



dqt qt 



3 |F(g)|2i?,(y)i?NB(gt 





— S(g,fi) / dv ^ 



g4|e(g)P 
v*{v + w) 



(CIO) 



Let us assume further that the dynamic structure factor S{q, Q) 
decreases rapidly with the growth of Q at low Q. Then we can 
set n = or w = everywhere in the integrand of (C7) or 
(CIO) except in S{q,Q). In this way we have Rt{jj,u,w) = 
RT{y,u) (cf. (20) and (C8)). The integration ove r Q. yields the 
static structure factor S{q) in accordance with (CI), and the 



Coulomb logarithm (C7) transforms into the static Coulomb 
logarithm (19). 

Thus Eq. (C7) extends the results of Sect. 3.3 to the case 
of the dynamic screening. In the latter case the Coulomb log- 
arithm is given by a three dimensional integral which could be 
computed if the dynamic structure factor were known. How- 
ever the detailed calculations of S{q, SI) have not yet been per- 
formed. We expect that the formalism will be useful in the 
future after S{q, Q,) is determined. Our preliminary estimates 
indicate that the dynamic effects are not very significant. 
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